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RESUMEN 

Presentamos un estudio astrometrico de la radiofuente VLA 1 de donde se origina el 
flujo HH 111 usando nuevos datos (2007) asf como tambien observaciones de archivo 
(1992-1996). Todos los datos fueron tornados a 3.6 cm con el Very Large Array en su 
configuracidn mas extendida (A). La fuente VLA 1 ha tenido un cambio dramatico, 
mostrando una eyeccion monopolar en la epoca de 2007. Tambien reportamos la 
deteccion de una fuente de continuo compacta (VLA 3) localizada a (-10'.'6, 98'.'7) 
de VLA 1 . No se encontraron movimientos propios significativos en VLA 1 ni en 
VLA 3 y los lmiites superiores son consistentes con los que se encontraron para 
radiofuentes (embebidas) en la Nebulosa de Orion. Favorecemos la interpretacion 
de que a 3.6 cm estamos observando dos chorros casi perpendiculares. HH 111 
representa un nuevo caso de eyecciones monopolares en un objeto joven estelar. No 
se pudo establecer si VLA 3 es Galactica (o extragalactica). 

ABSTRACT 

We present an astrometry study of the radio source VLA 1 at the core of the HH 111 
outflow using new data (2007) as well as archival observations (1992-1996). All data 
were taken at 3.6 cm with the Very Large Array in its most extended (A) configura- 
tion. The source VLA 1 has undergone a dramatic morphological change, showing 
a one-sided knot ejection in the 2007 epoch. We also report on the detection of a 
3.6 cm compact continuum source (VLA 3) located at (— 10'.'6, 98'.'7) from VLA 1. 
No significant absolute proper motions were found for VLA 1 and VLA 3 and the up- 
per limits are consistent with those found for (embedded) radio sources in the Orion 
Nebula. We favor the interpretation that in the continuum at 3.6 cm we are observing 
two nearly perpendicular jets. HH 111 presents a new case of one-sided jet ejection 
in a young stellar object. The Galactic (or extragalactic) nature of VLA 3 remains 
unclear. 

Key Words: Astrometry — ISM: individual objects: HH 111 — ISM: jets and out- 
flows — radio continuum: ISM 

1. INTRODUCTION 

Many observati ons have increasingly shown that stars are rarely born alone (e.g., 
Lada & Ladall2003l) . 3inary and multiple systems may be surrounded by circumstel- 



lar and/or circumbinary disks, and often drive episodic jets and powerful outflows. 
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For instance. iPech et al. | J2010|) presented Very Large Array (VLA) multi-epoch ob- 
servations of the very young hierarchical multiple system IRAS 16293-2422. These 
authors found a bipolar ejection and suggested it was a consequence of an episode of 
enhanced mass loss possibly produced by an increase in accretion onto the protostar. 

An interesting binary system is fo und at the core of the giant Herbig-Haro ob- 
ject HH 111, originally discovered bv lReipurthl(ll989|). This flow in Orio n has been 
extensively studied since its discovery (see iNoriega-Crespo et al. 201 li and refer- 
ences therein). It is highly collimated approximately in the east-west direction and 
displays a large number of individ ual knots moving in the plane of the sky with 
velocities of several hundred km s _1 dReipurth et al ll992h . We will refer to this east- 
west je t as the HH 111 optical jet. Near-infrared observations by Gredel & Reipurth] 
(119931) revealed the presence of a second outflow (HH 121) nearly perpendicular to 
the HH 111 one, making the system a quadrupolar outflow. We will refer to this 
second jet, approximately in the north-south direction, as the infrared jet. 

The driving sources of the HH 111 jets are believed to be the components of 
the suspected binary, class I sou rce IRAS 05491+024 7, which was associated to the 
radio source VLA 1 reported by iReipurth et al.l (1 1 999D at 3.6 cm. These authors also 
detected an additional faint, u nresolved, radio source (VLA 2) to the north-west of 
VLA 1. Recently. iLed (1201 lb imaged the HH 111 protostellar system at 1.3 mm. A 
faint (4cr) source appeared at the same position of VLA 2 possibly tracing its putative 
du sty disk. 



Rodriguez et alj (120081) presented observations of the system at 7 mm and dis- 



cussed three possible interpretations for VLA 1, since the emission revealed a struc- 
ture that can be described as two overlapping, almost perpendicular, elongated sources. 
They favor two interpretations in which the quadrupolar structure was either a disk 
elongated in the north-south direction with a perpendicular jet or two orthogonal 
disks. In the third interpretation, the emission comes from two jets. Nevertheless, 
the free-free contribution is only about 30%, hence the emission is dominated by 
dust. However, it is possible that the core is dominated by dust an d the extended 
emission dominated by free-free emission (see lRodrfguez et alJ l2008). The emission 
at 1.3 mm that showed a structure elongated in the nor th-south direction and perpen- 
dicular to the HH 1 1 1 optical je t ( Reipurth et aljfl999l) was attributed to the thermal 
emission from a disk jLeelboil ). 

In this paper we present sensitive, high angular-resolution 3.6 cm continuum 
observations at the core of the HH 111 outflow to search for absolute proper motions 
and to study the evolution of the HH 1 1 1 jet in VLA 1 . In Sect. [2] we describe our 
observations carried out with the Very Large Array as well as the data reduction. In 
Sect. [3] we present the results together with the analysis. The discussion is presented 
in Sect.|4]in the context of one-sided knot ejections versus bipolar ejections and the 
summary is given in Sect. [5] 



2. OBSERVATIONS AND DATA REDUCTION 



The 3.6 cm (v = 8.4 GHz) continuum observations of the core of the HH 111 out- 
flow were made with the Very Large Array (VLA) of the National Radio Astronomy 
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TABLE 1 

VERY LARGE ARRAY OBSERVATIONS AT 3.6 CM. 



Epoch 


Project 


Phase calibrator 


Bootstrapped flux density 


Synthesized beam" 


a 












(mJy beam ^ ) 


1992 Nov 2 (1992.84) 


AR278 


0541-056 


0.8691 ± 0.0002 


0'.'26x0'.'23;+42° 


0.02 


1992 Dec 18 (1992.96) 


AR278 


0541-056 


1.0253 ±0.0002 


0'.'25 x 0'.'24; +37° 


0.02 


1992 Dec 19 (1992.97) 


AR278 


0541-056 


1.0229 ± 0.0002 


0'.'28x0'.'25;+8° 


0.02 


1994 Apr 30 (1994.33) 


AR277 


0552+032 


0.7663 ± 0.0003 


0'.'21 x0'.'20;-57° 


0.02 


1996 Nov 11 (1996.86) 


AR367 


0552+032" 


0.8085 ± 0.0001 


0'.'25 x0'.'22;-12° 


0.01 


1996 Dec 28 (1996.99) 


AR277 


0552+032 


0.6997 ± 0.0001 


0'.'25 x 0'.'22; -5° 


0.02 


1996 Dec 29 (1997.00) 


AR277 


0552+032 


0.7088 ± 0.0001 


0'.'30x0'.'25;+3° 


0.02 


2007 Aug 12 (2007.61) 


AG747 


0552+032 


0.6780 ± 0.0003 


0'.'26x0'.'23;-ll° 


0.01 



"Major axis x minor axis; position angle of major axis. 
b For this epoch, the source 0541-056 was also observed. 



Observatory (NRAC0) in the A configuration on 2007 August 12. The phase center 
of the observations was at a(J2000) = 5 h 51 m 46?25, 5(J2000) = 2°48'29"6. The phase 
calibrator 0541-056 was observed for the first three 1992 epochs while the source 
0552+032 was used as phase calibrator for the rest of the epochs (see Table [TJ. The 
absolute amplitude calibrator was the source 1331+305. We used a model image for 
1331+305 provided by NRAO in order to improve flux accuracy. 

In addition, we used seven epochs taken from the VLA archive. The eight ana- 
lyzed data sets span about 15 years and are listed in TableQ] The data were analyzed 
in the standard manner using the Astronomical Image Processing System (AIPS) 
software package; the calibrated visibilities were imaged using weights intermediate 
between natural and uniform (with the ROBUST parameter set to 0) and CLEANed 
in an iterative fashion. To obtain accurate absolute astrometry, we precessed the data 
taken in B1950.0 to J2000.0 using the task uvfix and used the most recent position 
of the phase calibrator for all epochs. We expect that the residual systematic er- 
ror affecting the data to be around 10 milliarcseconds (mas). These residual errors 
were added in quadrature to the positional uncertainty delivered by a Gaussian fitting 
program (task jmfit of AIPS). 



3. RESULTS AND ANALYSIS 



We detected the source VLA 1 at epoch 2007.61 and show the image in the upper 
panel of Fig. [T] The morphology of VLA 1 changed dras tically between the early 
epochs feodriguez & Reipurthlll994l iReipurth et alJll999h and the 2007.61 epoch. 
In the last epoch there is a knot ejection to the east which was not present in previous 
observations. Parameters of the ejecta and discussion of the one-sided ejection are 
presented in Sects. I3~3l and l4. II respectively. Weak emission (at 3<x level) to the north 
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TABLE 2 
ABSOLUTE PROPER MOTIONS. 



Source 


Ha cos 5 


He 




(mas yr _1 ) 


(mas yr _1 ) 


VLA 1 


-0.34 + 0.50 


-0.75 + 0.93 


VLA 3 


-0.33 + 0.57 


2.41 +2.23 



and south of VLA 1 is also detected. We believe this weak emission is associated 
with the base of the infrared jet. 

In the 2007 observations, due to potential problems with several antennas that 
were in transition to the K. Jansky Very Large Array (VLA) system, the rms noise 



level w as not sufficient to detect the VLA 2 source, as previously reported in Reipurth et al 



(Il999h . 

On the other hand, we report on the detection of a radio source (hereafter VLA 3) 
located at (-10'.'6, 98'.'7) from VLA 1. Inspection of the maps at other epochs con- 
firms the presence of this source. After correcting the image for the primary beam 
response, VLA 3 has a flux density of S — 0.22 ± 0.03 mJy in the 2007 image. The 
flux density at all four epochs is consistent with a value of S = 0.3 + 0.1 mJy. There 
is a suggestion of time variability, but the data of 1992 and 1994 do not have a suffi- 
cient signal-to-noise ratio to reach a firm conclusion. The elongation of the source in 
the north-south direction seen in the bottom panel of Fig. Q] is caused by bandwidth 
smearing and is not intrinsic to the source. 

3.1. Proper motions 

In order to search for absolute proper motions of VLA 1 and VLA 3, we con- 
catenated the epochs 1992.84, 1992.96, and 1992.97 into epoch 1992.92; and the 
epochs 1996.86, 1996.99, and 1997.00 into epoch 1996.95. Figure|2]shows the radio 
positions of VLA 1 and VLA 3 as a function of time. We performed a least squares 
fit to four epochs (1992.92, 1994.33, 1996.95, and 2007.61) and the results are listed 
in Table[2] We found that VLA 1 and VLA 3 do not show significant absolute proper 
motions and they are consist ent within 2cr with those found for (embedded) radio 
sources in the Orion Nebula dGomez et al 1 l2005h . 



Further observations will be necessary to detect VLA 2, so that we are able to 
obtain relative proper motions with respect to VLA 1 and test the hypothesis that 



these sources form a non-hierarchical triple system in disintegration dReipurfh et al 



1999). 



3.2. On the nature of VLA 3 

In order to obtain the expected number of extragalactic sources, N, in the field 
of study with a flux d ensity comparable or larger than that of VLA 3, we follow 



Windhorst et all d 19931) 



-3 M P 



N = 2.4 x 10^ arcmin" 2 , (1) 

\mJy/ 



A ONE-SIDED KNOT EJECTION IN HH 1 1 1 



5 



VLA 



S 2 o 48'30" 



2°48'29 




5 h 51 m 46 s 3 



tb.2 



Right Ascension (J2000) 



VLA 3 



2°50'09 



124 AU 



2 D 50'08 M 




to phase center 



5 51 4576 45.5 
Right Ascension (J2000) 



Fig. 1. Contour images of the 3.6 cm continuum emission of VLA 1 {top) and VLA 3 {bottom) 
at epoch 2007 .6 1 . The synthesized beam (0'.'26 x 0723 ;PA= -11°) is shown in the bottom left 
corner. Top: First contour and contour spacing for VLA 1 are 31.5 /jjy beairr 1 (3cr). The ar- 



rows in the east- west direction indicate the P. A. of the HH 1 1 1 jet (97.5° uNoriega-Crespo et al 



201 lh . The arrows in the north-south direction indicate the orientation of the HH 121 infrared 
jet. Bottom: First contour and contour spacing for VLA 3 are 42.0 yuJy beam 1 (3<r). The bar 
indicates the linear scale. The arrow points in the direction to the phase center. The map has 
been corrected for the primary beam response. 
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Fig. 2. Absolute proper motions in right ascension (or) and in declination (S) for VLA 1 {upper 
panels) and the newly reported radio source VLA 3 (bottom panels) as a function of time. The 
dashed line in each panel represents the least squares fit to the data. The right ascension axis 
is given in seconds, while the declination axis is given in arcseconds. 
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where S is the flux density in mJy. For S - 0.3 mJy and a field of ~ 3.4x3.4 arcmin 2 , 
N corresponds to 0.13 sources. Thus, VLA 3 only has a 13% probability of being 
extragalactic. 

VLA 3 does not show significant proper motions, as expected for an extragalactic 
source, but is also consistent with Orion sources. Therefore, the measurement of 
proper motions, in this c ase, is not useful in determining its nature. Very recently, the 



AMI Consortium et al.l (12012b detected a radio source whose emission peaks at the 
position of VLA 3 at 1 .9 cm (16 GHz) with a beam size of 69" x 27" . They obtained 
a spectral index, a ami, of -0.23 + 1.10, a value consistent with those found for 
extragalactic sources, but whose large error makes it consistent with several different 
emission mechanisms. 

We s earched for possible c ounterpart s in three point source catalogs, na mely, the 



AKA RI dMurakamietalJ2007l) . 2MASS dskrutskie et alkood) . and WISE dwrirfit et al 



201 Oh catalogs. We found no association with sources from the AKARI and 2MASS 



catalogs. A WISE point source was found as possible counterpart of VLA 3 at a 
distance of The source is detected in the 3.4 fim and 4.6 //m bands. There are 
marginal signals in the 12 //m and 22 //m bands that are considered as spurious or 
affected by other artifacts in the catalog. 

Although VLA 3 is more likely an extragalactic source, in the remaining of this 
contribution we restrain ourselves from further interpretation of VLA 3. 



3.3. VLA 1: ejecta parameters 



The parameters of the ejecta are estimated f rom observables at epoch 2007.61. 
Using the model by Mezger & Henderson! (119671) . and assuming optically thin free- 
free emission, the mass of ionized gas, M,, and the electron density, n e , of a homo- 
geneous sphere can be calculated as 



Mi 



I T \ 0175 / v \0 05 I c y 
3.39 x 10-5 (^—) — 

\10 4 K/ UGHz/ \lmJy/ 

, T ,0.175 , >0.05 / C \°- ; 

1 kpc / \ 1 " / 



(2) 



(3) 



where T e is the electron temperature, v is the frequency, S v is the observed integrated 
flux density, D is the distance, and 6 is t he angular size. We assume a T e of 10 4 K 
at 8.4 GHz, and use a distance of 414 pc ( Menten et alj|2007 ). The flux density was 
obtained by fitting a Gaussian to the ejecta. The emission was resolved in only one 
direction with a deconvolved angular size of 0'.'28, while in the other direction the 
emission comes from a region smaller than 0'.'23. Using the fitting parameters of 
S =0.18 mJy and the geometric mean of the deconvolved angular size of < 0'.'25, 
we obtain M, < 2.2 x 10~ 7 M and n e > 4.2 x 10 4 cm" 3 . 
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4. DISCUSSION 



4.1. A one-sided knot ejection 

Figure [3] compares the images for VLA 1 of four epochs, including the 2007 
observations. We show images of concatenated data for ep ochs 1992 and 1996 (see 
Sect. 13. 11 1. For the 1992 epoch. iRodriguez & Reipurthld 19941) note that the emission 
of VLA 1 consists of a single component elongated to the west, possibly tracing 
a faint ejection. The elongation is to the east in the 1994 epoch, and to the east- 
west direction for the structure seen in the 1996 epoch, whereas in the new 2007 
observations, we additionally see a second component to the east, which we interpret 
as a one-sided knot ejection (see also Fig. [TV Moreover, weak emiss ion is detected 
in the north-south direction, and according to Rodriguez et al . (2008), it probably is 
du e to the presence of th e second (infrared) HH 121 jet. In the 3.6 cm map presented 
bv lReipurth et al.l(11999l) . VLA 1 shows clear elongations in the north-south direction 
after using a maximum entropy reconstruction (see their Fig. 3). 

It is important to point out that the map corresponding to the 1996 epoch, which 
has a better rms compared to the 2007 epoch, presents an elongation to the east but 
that the ejection cannot be resolved. It is also inter esting that, in the o ptical, the 
(eastern) redshifted lobe is completely obscured (see lReipurth et al.ll 19921) . contrary 
to what we see at radio wavelengths in the 2007 epoch where the ejection is clearly 
seen to the east. 

Dust extinction affects observations at optical and infrared wavelengths, while 
the radio regime suffers the le ast. One may thi nk that asymmetries often seen in 
jets at the former wavelengths (IHirth et al.lll994l) are mainly due to the material in 
the parental cloud, where the protostars are embedded, while observations in the 
radio would show symmetric jets. In reality, this is not alwa ys the case, for ins tance, 



Curiel et al. 



(2006) see 



for the (high-mass) young stellar object Cepheus A HW2 
an asymmetry in the radio lobes and explain it by means of proper motions of the 
exciting source and by temporal variations in the velocity of the ejections. As we 
mentioned in Sect. 13.11 VLA 1 does not show significant absolute proper motions, 
hence the scenario in which the driving source has proper motions can be ruled out. 
New and more sensitive VLA observations will be necessary to also measure proper 
motions of the ejecta. An other possibility for the ejection is due to episodic mass 
loss rate increments (e. g.. IPech et al. 201ob. Ot her examples of asymmetry in the 
radio were presented bvlRodrfgue z et al.1 (l2012alfbT) for the (low-mass) stars DG Tau 
and DG Tau B, where similar one-sided knot ejections were noticeable in the 3.6 cm 
maps. 

There are two model s for describing the launching region of jets: the disk- winds 
(IKonigl & Pudrita 12 000) when they originate from the accretion disk and the X- 
winds ( Shu et alJuOOOl) when they originate in the disk-star in terface. Although a 
lot of progress has been done in star-disk simulations (see, e.g., Pu dritz et al.ll2007l) . 
questions of whether these winds co-exist or one is at work while the other is not 
still remain unanswered. Possible causes for jet asymmetries may be traced back to 
asymmetries in the disk, e.g., a warped disk. Obviously, further theoretical work on 
the mechanisms at work in episodic, one-sided jets is urgently needed. 
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Fig. 3. Contour images of the 3.6 cm continuum emission of VLA 1 at different epochs. 
Details for epoch 2007.61 are the same as in Fig. [T] First contour and contour spacing 
for VLA 1 are, from top to bottom, 36.0 //Jybeairr 1 (3cr), 54.3 yijybeairr 1 (3cr), and 
21.0 yuJybeairT 1 (3cr), for epoch 1992.92, 1994.33, and 1996.95, respectively. The synthe- 
sized beams are shown as filled ellipses for each epoch [(0'.'26x0'.'24; PA = 34°), (0'.'22x0'.'21; 
PA = -41°), (0'.'25 x 0'.'22; PA = -8°)]. 
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Observations of VLA 1 in HH 1 1 1 together with those of DG Tau, DG Tau B, 
and Cepheus A HW2 suggest one-sided knot ejection events do occur in young stellar 
objects. 

5. SUMMARY 

We have performed 3.6 cm continuum observations with the VLA together with 
the analysis of additional archival data toward the core of the HH 111 outflow. Our 
main results can be summarized as follows: 

• The source VLA 1 has undergone a dramatic morphological change, showing 
a one-sided knot ejection at epoch 2007. 

• We have confirmed that the base of the infrared jet HH 121 is detectable at 
centimeter wavelengths, although it is weaker than the emission from the base 
of the optical HH 1 1 1 jet. 

• We have reported on the detection of the radio source VLA 3 located to the 
north of VLA 1 . However, its Galactic or extragalactic nature remains unclear. 

• We have found that VLA 1 and VLA 3 do not show significant absolute proper 
motions and the upper limits determined by us are consistent with the proper 
motions found for (embedded) radio sources in the Orion Nebula. 

L. G. was supported for this research through a stipend from the International 
Max Planck Research School (IMPRS) for Astronomy and Astrophysics at the Uni- 
versities of Bonn and Cologne. L. L. and L. F. R. acknowledge the financial support 
of DGAPA, UNAM and CONACyT, Mexico. L.L. is indebted to the Alexander von 
Humboldt Stiftung and the Guggenheim Memorial Foundation for financial support. 
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